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We investigate J/ip transverse-momentum (pt) distributions and their centrality dependence in 
heavy-ion collisions at SPS and RHIC within the framework of a two-component model, which 
includes (i) primordial production coupled with various phases of dissociation, (ii) statistical coa- 
lescence of c and c quarks at the hadronization transition. The suppression of the direct component 
(i) is calculated by solving a transport equation for J/ip, Xc and ip' in an expanding fireball using 
momentum dependent dissociation rates in the Quark-Gluon Plasma (QGP). The coalescence com- 
ponent is inferred from a kinetic rate equation with a momentum dependence following from a blast 
wave approach. At SPS energies, where the direct component dominates, the interplay of Cronin 
effect and QGP suppression results in fair agreement with NA50 pt spectra. At RHIC energies, the 
Pt spectra in central Au+Au collisions are characterized by a transition from regeneration at low p t 
to direct production above. At lower centralities, the latter dominates at all pt- 

PACS numbers: 25.75.-q, 12.38.Mh, 14.40.Lb 
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I. INTRODUCTION 

It has long been suggested that the suppression [l[ of 
J/ip mesons can be utilized as a probe of Quark-Gluon 
Plasma (QGP) formation in ultrarelativistic heavy-ion 
collisions (URHICs). This effect has indeed been ob- 
served in Pb-Pb collisions at the Super Proton Syn- 
chrotron (SPS) [2j, consistent with theoretical models in- 
voking charmonium dissociation reactions in a QGP 
H®. At collider energies, however, a copious production 
of cc pairs has led to the suggestion that charmonia can 
be regenerated by a coalescence of c and c quarks close 
to the hadronization transition [H, 0, [!]• The presence 
of this mechanism is supported by data from the Rela- 
tivistic Heavy-Ion Collider (RHIC) Q, where, despite the 
higher temperatures of the putative QGP, the observed 
suppression is similar to SPS energies (as, e.g., predicted 
in Ref. @). However, a quantitative assessment of the 
regeneration and suppression mechanisms at RHIC has 
not been achieved yet. 

Transverse momentum (p t ) spectra of charmonia are 
hoped to provide additional means of discrimination. In 
Ref. the average (pf) of primordial J/ip's with QGP 
suppression has been computed, while in Refs. [ljj [H| 
the regeneration component has been studied within a 
blast wave description based on thermalized J/ip mesons 
at the QCD phase boundary. In Refs. [H,[l3| the impact 
of non-thermalized c quark distributions on the pt spec- 
tra resulting from recombination has been studied, but 
no direct component was accounted for. In Ref. [14( the 
average (pf) and elliptic flow at RHIC has been calcu- 
lated including both direct and regenerated components 
using gluo-dissociation rates with vacuum binding ener- 
gies for the J/ip and Xc- 

In the present Letter we provide a comprehensive de- 
scription of p t spectra at both SPS and RHIC, includ- 
ing their centrality dependence and absolute yields. A 
proper description of the inclusive yields vs. centrality 



turns out to be particularly important for interpreting p t 
spectra in terms of direct and regeneration contributions. 
We adopt a previously constructed two-component ap- 
proach 3] employing inelastic charmonium reaction rates 
(extended to finite 3-momentum) which account for re- 
duced binding energies in the QGP. The latter point is 
essential to allow for a realistic treatment of \c and tp' 
states, which are expected to be measured in the future 
but also make up 30-40% of the inclusive initial J/ip yield. 
In the following, we recall the basic ingredients of the 2- 
component model and its extension to finite momentum 
(Section[TT| , apply our approach to heavy-ion collisions at 
SPS and RHIC (Section EH) and conclude (Section [TV} . 



II. TWO-COMPONENT APPROACH 

In analogy to the case of light hadrons, one may de- 
compose pt spectra of charmonia (\P = J/ip,Xo VO m 
URHICs according to their production mechanism into 
hard (high-p t ) and soft (low-pt) components, 
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where the direct component (first term) is associated 
with hard production in primordial N-N collisions, sub- 
ject to suppression in the subsequent medium evolution. 
The soft component (second term) is conceptually sim- 
pler than in the light sector, based on the notion that 
c and c quarks are exclusively produced primordially, 
leaving their coalescence as the only source of secondary 
charmonium formation. Since regeneration is governed 
by the phase space density of c and c quarks in the 
medium, it is sensitive to the c-quark momentum spec- 
tra. Indirect measurements of c- quark sp ectra at RHIC 
(via semileptonic decay electrons) [15l . Ilq ] indicate strong 
rescattering effects which in theoretical models [17[ imply 
an approximate thermalization up to c-quark momenta 
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of pi ~ 2 — 2.5 GeV [lj|. Thus, following our earlier 
developed 2-component model [H, [l!|, we approximate 
the coalescence component with a thermal blast wave 
description, while direct production is computed in a mi- 
croscopic suppression calculation in QGP and hadronic 
phase, as will be detailed in the remainder of this sec- 
tion. Both terms in Eq. JT]) are evaluated in the same 
expanding fireball model. 

Let us first address the direct component. Since 
the charmonium masses are much larger than the typ- 
ical temperature of the medium, a Boltzmann transport 
equation is appropriate to describe the time evolution 
of the phase space distribution, f\&(x,p, r), through the 
QGP, mixed and hadron gas (HG) phase, 

P^d,j,fa(x,p,T) = -Ey r^(x,p,T) fa(x,p,r) , (2) 



with 



where = y/ to|, + p 2 is the energy of \& with 3- 
momentum modulus p, and x is its position in the fireball. 
For simplicity, we constrain our calculation to the longi- 
tudinal rest frame of the charmonium 0], i.e., solve the 
Boltzmann equation in (2+l)-dimensions. For the ini- 
tial charmonium distribution we assume a factorization 
into spatial and momentum dependencies, f(xt,pt, To) — 
f{xt,To)-f(pt,To), where To is the initial (thermalization) 
time of the medium (QGP or mixed phase). The spatial 
part of the initial distribution is obtained from a Glauber 
model including nuclear absorption, 

fa(x t ,T ) = cr* / d 2 s dz dz'p A (s,z) p B (x t - s,z') 



dz A p A (s, Z A )<Jr, 



x exp 



dz B PB(x t - S,Z B )cTnuc > , (3) 



x exp 



where p A ,B ar e Woods-Saxon profiles [2(| of nuclei A 
and B and cr* is the VP production cross section in p+p 

collisions (we use dap^/dyiy = 0)^25(750) nb at SPS 
(RHIC) [U HH). The nuclear absorption cross sec- 
tion, cr nuc , serves as a parameter to account for pre- 
equilibrium charmonium suppression due to primordial 
nucleons passing by, estimated fro m p -A collisions; at 
SPS we adopt the values of Refs. [H,[24j], <j nuc =AA mb for 
J/ip, Xc and 7.9 mb for tp' . At RHIC, we use a nuc =1.5 mb 
(for J/ip, Xc) based on Ref. [25[, which is compatible 
with a recent update [26| (for a nuc =2.7 mb the total J/tp 
yield in our model decreases by 8% for central Au-Au at 
RHIC). For -0', we employ an accordingly increased value 
of 2.7 mb. The initial momentum spectra are obtained 
from p+p data, augmented by a Gaussian smearing to 
simulate nuclear p t -broadening (Cronin effect), 



fa (ft, to) 



27TCT 2 



d 2 p' t exp(--^) f N N(\pl -Pt\) 



(4) 

where fjsrNipt) is the spectrum in elementary N-N 
collisions. At SPS, f NN (Pt) = exp(-pf /(pf)) 



(p 2 t ) = 1.1 5 GeV 2 /c 2 [13, H|, and at RHIC 
fN N (pt) = A (1 +p 2 /B 2 )- 6 with 5=4.1 GeV 2 yielding 
(p|) P j,=4.14 GeV 2 [22j |. The Cronin effect is computed 



/pa- 
using 2a 2 =a g N ■ (I) where (I) represents the centrality 



dependent mean nuclear path length of the gluons be- 
fore fusing into \E' [2^]. At SPS, the extracted coefficient 
is a s jv=0.076GeV 2 /fm (HHH, while a fit to d-Au data 
at RHIC gives a g N— 0.1 GeV 2 /fm with a rather large 
uncertainty (e.g., a s jv=0.6 GeV 2 /fm is still compatible 
with d+Au data, but results in R AA (p t =6 GeV)~9 be- 
fore QGP suppression in 0-20% central Au+Au). 

The most important microscopic ingredient to the 
transport Eq. ^ are the charmonium dissociation rates, 
r*, which can be expressed via inelastic cross sections, 
a^f s , for ^ scattering on medium constituents i as 



r*(z, p, r) J 



d 3 k 
(2^)3 



/*(w fc ;T(r)) <r v rel (5) 



with v rei = F/(EyEi), E t = (k 2 + m 2 ) 1 / 2 , flux fac- 
tor F — ((p^kfj,) 2 — m^rnfy 1 / 2 , k^: parton/meson 4- 
momentum, f l (ujk',T): thermal Fermi/Bose distribution. 
In the QGP, color Debye screening is expected to re- 
duce charmonium binding energies, e B (which eventu- 
ally vanish), which is supported by recent lattice QCD 
calculations [3(J [H|. Under these circumstances, gluo- 
dissociation reactions, \& + g — > c + c, become ineffi- 
cient and should be replaced by quasifree dissociation, 
i + 'fy^i + c + c (i=g,q,q) Q . Here we extend these cal- 
culations to finite 3-momentum and compare the rates 
to gluo-dissociation (using vacuum binding energies and 
vanishing thermal gluon mass) [32] in Fig. [TJ The strong 
coupling constant a s in the quasifree cross section is 
one of two adjustable parameters in our approach and 
is fixed to reproduce the J/ ip yield in central Pb-Pb col- 
lisions at the SPS (resulting in a s =0.24). Except for 
J/ip's at rather low temperatures and 3-momenta, the 
gluo-dissociation rate decreases with increasing p due to 
a pronounced maximum in the pertinent cross section 
(at a gluon energy to ~ 1.43e^ in the rest system of 
the J/ip [111]). On the other hand, the quasifree rate al- 
ways increases with p due to a smoothly increasing cross 
section with center-of-mass energy, similar to Ref. [34j . 
This reiterates the importance of using the quasifree rate 
(rather than gluo-dissociation) for small binding ener- 
gies, especially at finite momentum. The increase of the 
quasifree rate with p is more pronounced at low tempera- 
ture since at high temperature most partons are energetic 
enough to destroy a J/tp irrespective of its momentum. 
This trend is weaker for the Xc (lower panel) due to its 
small binding energy: even at low temperature most par- 
tons carry sufficient energy to destroy it. 

In the HG, we employ inelastic cross sections with ir 
and p mesons from a flavor-S'?7(4) effective Lagrangian 
approach (35l.[36j. 

The final ingredient required to solve the transport 
equation @ is the space-time and temperature evolu- 
tion of the system, which we model by an isentropically 
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FIG. 1: (Color online) Comparison of the momentum de- 
pendence of quasifree (solid lines) and gluo-dissociation rates 
(dashed lines) for J/ip (upper panel) and \ c (lower panel) at 
different temperatures. 



expanding fireball model represented by a cylindrical vol- 
ume, 

Vfb(t) = {z + v z t + ^a 2 r 2 ) tt (r + ^aj_T 2 ) 2 , (6) 

where Zq, v z , a z , tq and a± are the initial longitudinal 
length, longitudinal expansion velocity and acceleration, 
initial transverse radius and transverse acceleration, re- 
spectively. At fixed total entropy (matched to the ob- 
served hadron multiplicity at given centrality and colli- 
sion energy), the temperature of the system follows from 
the equation of state (massive partons in QGP and reso- 
nance gas in HG) . We update the transverse acceleration 
to a ± =0.08 c 2 /fm and 0.1 c 2 /fm at SPS and RHIC as 
used in recent applications of the fireball model to dilep- 
ton [32j and heavy-quark [l7| observables. The initial 
temperature for central A- A collisions at SPS (RHIC) is 
210(370) MeV, with thermal freezeout at T {o ~ 110 MeV. 

The Boltzmann transport equation ^ can now be 
solved for the final phase-space distribution f^(x t ,pt, t/) 
of J/ip, Xc and ip 1 at the freeze-out time ry for fixed cen- 
trality and collision energy. Upon integration over the 
transverse plane of the medium we obtain the pt spec- 
trum of the direct component as 



dA* 



Ptdpt 



= Id x t f(x t , pt, Tf) , 



(7) 



It is worth noting that the leakage effect [38| is imple- 
mented by switching off the suppression if a charmonium 
state moves outside the fireball, i.e., r*(af, p, t) =0 for 
Xt{r) > ro + ^a±r 2 . As we will see below this effect is 
significant for charmonia at high p t . 

Let us now turn to the coalescence component, i.e., the 
second term on the right-hand side of Eq. {1} . As stated 
above, we assume the regenerated charmonia to follow 
a local thermal equilibrium distribution with transverse 
flow velocity given by the blastwave expression |39j . 



Ptdpt 



ex m t 
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m t coshyt 
T 



p t sinh y t 
T 



(8) 

(mt=y/ rriy + p 2 ). Since charmonium regeneration is in- 
operative in the HG [H, E3] > we evaluate the blast wave 
formula at the hadronization transition as following from 
the fireball model, Eq. ©, with T=T C =170(180) MeV 
and transverse rapidity y t =tanh~ 1 v t (r) using a lin- 
ear flow profile Vt(r)=v s j^ with surface velocity 
v s =0.33(0.49)c and transverse radius i?,=7.3(7.9) fm for 
central collisions at SPS (RHIC). To determine the nor- 
malization of the coalescence component we utilize a 
momentum- independent rate equation [4fJ], 



diV* 

~d~r 



(9) 



where = T>i>(p = 0) and is the equilibrium num- 
ber of charmonia for a given number of cc pairs in the 
system (based on total cross sections (7^=5.5(570) fib 
at SPS (RHIC)). The charmonium yield due to the gain 
term is identified with the abundance of the coalescence 
component. As in Refs. [lj| H(| a thermal relaxation 
time, m , for charm quarks is introduced to mimic 
a reduced charmonium equilibrium limit due to incom- 
plete kinetic equilibration via a relaxation factor 71=1- 
exp(— J G?r/r* hcrm ). Due to the current uncertainties in 
<jpp, r' herm and its schematic implementation we cannot 
quantitatively predict the coalescence yields. Therefore 
we adjust r* herm to the inclusive J/ip yield in central 
Au-Au at RHIC. This point will be improved in future 
work by solving the rate equation at finite p with (time- 
dependent) c-quark momentum distributions as obtained 
from Langevin simulations [TtJ which result in fair agree- 
ment with the semileptonic single-electron Raa and x>2 at 
RHIC [H [3. Here, we employ r c therm =7 fm/c (in line 
with the microscopic approach of Ref. [4l|). compared to 
^2-4 fm/c in Ref. [l/J. This update reduces the regen- 
eration yield by 30-50% in central collisions at SPS and 
RHIC (see Sec. HI] below). 



III. J/ip YIELDS AND SPECTRA AT SPS AND 
RHIC 



We start our phenomenological analysis at SPS ener- 
gies. The centrality dependence of inclusive J/ip produc- 
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FIG. 2: (Color online). NA50 data @] for the central- 
ity dependence of J/i/;/Drell-Yan dimuons at SPS com- 
pared to our results with (black line) and without (purple 
line) 3-momentum dependence (both lines essentially coin- 
cide). The sensitivity of the coalescence component to the 
charm-quark relaxation time is indicated by the dash-dotted 
(r c thorm =7 fm/c) and dash-double-dotted (r* herm =3 fm/c) 
lines. 
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FIG. 3: (Color online), (pf) as a function of centrality at SPS. 
NA50 data [2?], Hs| are compared to our model calculations. 
The {pi) for the direct component (dashed line) and coales- 
cence component (dash-dotted line) are compared to the (p^) 
with nuclear absorption only including Cronin pt-broadening 
effect (dotted line). 



tion (including feeddown) is summarized in Fig. O To re- 
cover previous results from the momentum-independent 
calculations for central collisions, a minor reduction of 
a s from 0.26 to 0.24 in the quasifree dissociation rate 
has been applied (since the rate increases with p) . After 
this adjustment, there is no visible modification in the 
inclusive centralit y d ependence left compared to the pre- 
vious p=0 results [4(| (an increase of the leakage of J/ip's 
for smaller system sizes is essentially compensated by a 
decrease in the fireball lifetime). The sensitivity to the 
coalescence contribution is small, but this will be differ- 
ent at RHIC. Our calculated J/4> Pt spectra are used to 
compute its average {pf ) as a function of centrality, and 
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FIG. 4: (Color online). Results of the 2-component model 
for R-^{N part ) at RHIC, compared to PHENIX data [|. 
The previous (momentum-independent) model [40[, corre- 
sponds to the dash-double-dotted line (r* herm =3 fm/c) and 



the double-dash-dotted (t c 



=7 fm/c) line. The latter 



closely coincides with the total (solid line) using p-dependent 
rates (also with r* herm =7 fm/c). The dotted line represents 
the primordial input, while the dashed line additionally in- 
cludes QGP and HG suppression. The coalescence yield for 
-rthcrm =7 fm/c is given by the dash-dotted line. 



compared to NA50 data [23, [23 in Fig. El Most of the 
observed pt dependence follows from the Cronin effect 
of the primordial component, represented by the dotted 
line. The QGP suppression, which is stronger at high p t 
due to the increase of the dissociation rate with p (re- 
call Fig. [I]), leads to a slight reduction of (p^), improving 
the agreement with data. The coalescence component is 
rather insignificant. 

Next, we proceed to the centrality dependence of inclu- 
sive J/ip production in Au-Au collisions at RHIC, as rep- 
resented by the nuclear modification factor, RAA{N par t) 
(the number of J/^'s for a given number of participant 
nucleons, N partl relative to that in p-p collisions mul- 
tiplied by the number of binary collisions), cf. Fig. 01 
Previous p=0 results with the updated nuclear ab- 
sorption cross section (but with identical coalescence con- 
tribution) overestimate the most recent PHENIX data 
for A?p art >200; increasing r' hcrm to 7 fm/c improves this 
part at the expense of more peripheral collisions. The in- 
clusion of the 3-momentum dependence (with a s =0.24) 
does not resolve this potential discrepancy, despite the 
presence of the leakage effect, for similar reasons as de- 
scribed above for SPS energies. We note that the roughly 
equal partition of the 2 components for central collisions 
is quite similar to the results of Ref. [3| (where the vac- 
uum gluon dissociation mechanism has been employed). 

The key point is now how the inclusive centrality de- 
pendence of the two components reflects itself in the pt 
spectra. Our results for Raa{pi) for different centrality 
selections (approximated by the average number of bi- 
nary N-N collisions) is compared to PHENIX data 
in Fig. [51 As anticipated, the coalescence contribution 
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FIG. 5: (Color online). Raa vs. transverse momentum for 
different centrality selections of Au-Au at RHIC. PHENIX 
data [t| are compared to our model calculations: initial pri- 
mordial component (dotted line), including QGP and HG 
suppression (dashed line), and with leakage effect switched 
off (dash-double-dotted line); the coalescence contribution is 
given by the dash-dotted line. 



is concentrated at low pt (up to 2-4 GeV), most notably 
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FIG. 6: (Color online). Centrality dependence of (pf) within 
the 2-component model at RHIC, compared to PHENIX 
data 0]. The dotted line corresponds to the primordial in- 
put distribution (with nuclear absorption and Cronin effect), 
the dashed line includes QGP and HG suppression, and the 
dot-dashed line represents the coalescence component. 



for central collisions but quickly ceasing for more periph- 
eral ones. For the direct component, the pt-broadening of 
the Cronin effect induces an appreciable rise of RAA{pt) 
in the region from 2-5 GeV (cf. the dotted lines). This 
trend is largely counter-balanced by the QGP suppres- 
sion (dash-double-dotted line, with leakage effect off, im- 
plemented by ignoring the spatial fireball boundary) , due 
to the increase of the quasifree dissociation rate with mo- 
mentum. The leakage effect finally restores significant 
strength at higher p t (up to ~40% for p t >5 GeV) in the 
direct component (dashed lines) . The opposite p t depen- 
dence of the direct and coalescence spectrum combines 
into a rather flat total R,AA{pt) which is quite compatible 
with experiment. We emphasize that a proper descrip- 
tion of the absolute yields is an important ingredient to 
this finding (the underestimate for 20-40% central colli- 
sions could be improved upon with a somewhat smaller 
c-quark thermalization time). E.g., a pure coalescence 
spectrum can be compatible with the central data, but 
it would be less convincing for more peripheral collisions. 
Therefore, R.AA(pt) data may indeed discriminate a two- 
component from a one-component model, especially if the 
experimental uncertainty can be reduced. 

We have checked that within the current experimen- 
tal accuracy of Raa{pi) it is not possible to exclude 
different suppression mechanisms in the QGP medium. 
The data are also consistent with calculations employing 
a dissociation rate based on gluo-dissociation [32| with 
vacuum binding energy and zero thermal gluon mass, 
since the relevant suppression regime is for temperatures 
T < 300 MeV where the J/ i/> rate is only weakly momen- 
tum dependent. 

We finally condense the p t spectra into a centrality de- 
pendence of (p^), as computed from our spectra and com- 
pared to PHENIX data in Fig. [5] This plot reiterates 
the importance of the soft coalescence spectra in central 
collisions to provide a near-flat centrality dependence. 
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We recall, however, the currently large uncertainty in 
the Cronin effect as inferred from d-Au collisions. 



IV. SUMMARY AND CONCLUSIONS 

We have studied the 3-momentum dependence of J/ij) 
production in heavy-ion collisions based on a previously 
developed two-component model which accounts for pri- 
mordial production with subsequent suppression and sec- 
ondary regeneration close to the QCD phase boundary. 
For the direct component, we adopted a transport ap- 
proach including up-to-date empirical input for nuclear 
absorption and a Cronin effect in the initial state. The 
key microscopic ingredient is the charmonium dissocia- 
tion rate in the QGP. We argued that the quasifree de- 
struction mechanism provides a realistic treatment for 
small (in-medium) binding energies and the extension to 
finite 3-momentum. For the coalescence component, we 
adopted a blast wave description at the hadronization 
transition within the same fireball model used for the 
direct spectra. Our approach has essentially two param- 
eters: the strong coupling constant in the quasifree rate, 
which we adjust to the suppression in central Pb-Pb at 
SPS, and the thermal relaxation time of c quarks, which 
controls the magnitude of the coalescence component, 
adjusted to the J/tp yield in central Au-Au at RHIC. 
Within reasonable values for these parameters, a s ~0.24 
and r* herm ~5-7 fm/c, an approximate overall description 



of the centrality dependence of inclusive J/tp production 
at SPS and RHIC emerges. The key point of our Let- 
ter is that, without further assumptions, the calculated 
Pt spectra are largely consistent with available SPS and 
RHIC data. We argued that this supports the under- 
lying momentum dependence of the dissociation rate in 
connection with reduced binding energies, as well as the 
presence of a ^50% coalescence contribution in central 
Au-Au collisions at RHIC. More work is required to scru- 
tinize these findings, e.g., an extension to NA60 data at 
SPS, forward rapidities at RHIC, predictions for LHC, 
as well as more accurate input from d/p-A experiments. 
A microscopic transport treatment with c-quark spectra 
constrained by open-charm observables will be pursued 
and used to predict elliptic flow. Ultimately, the un- 
derlying charmonium properties in the QGP should be 
consistent with lattice QCD results to establish model- 
independent connections between the QCD phase dia- 
gram and the matter created in heavy-ion collisions. 
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